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Reliable predictions of cross sections and final-state distributions for 
QCD processes are a crucial ingredient in high-energy collider experiments, 
not only as a test of QCD but also for new particle searches. All system¬ 
atic approaches to this problem are based on fixed-order (FO) results in 
perturbation theory, and yield (usually at the next-to-leading order, NLO) 
the best available results for sufficiently inclusive observables. However, in 
many cases a more exclusive description of final states is needed. In such 
cases, in which one also combines perturbative calculations with a model for 
the conversion of partonic final states into hadrons, Monte Carlo (MC) sim¬ 
ulations are generally adopted. MC’s operate on partonic states with high 
multiplicity and low relative transverse momenta, which are obtained from a 
parton shower or dipole cascade approximation to QCD dynamics. This has 
to be confronted with FO results, which can describe the complementary 
region of small multiplicities, and large relative transverse momenta. 

The lack of large transverse momentum emissions, and the fact that total 
rates are computed to leading order accuracy only, are serious problems in 
MC simulations, especially when the CM energies are in the TeV range. 
These problems can be solved by a suitable combination of MC and FO 
methods. Given the flexibility of MC’s, it is actually desirable to use some 
FO techniques in the framework of MC simulations. One approach (matrix 
element corrections) acts at the level of the generation of hard processes 
in the MC: 2 —> n processes, with n > 2, are given to the parton shower 
as initial conditions (in standard MC’s, hard processes have usually a2-> 
2 kinematics). ME corrections therefore improve the description of large 
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transverse momentum emissions, but results for total rates have still a LO 
accuracy. 

A more involved approach to the problem aims at improving also the 
computation of total rates. There is a considerable freedom in the very def¬ 
inition of such an approach. Here I shall follow the formalism presented in 
ref. |], denoted by MC@NLO there. A MC@NLO is required to have the 
following features: a) The output is a set of events, which are fully exclu¬ 
sive. b) Total rates are accurate to NLO. c) NLO results for all observables 
are recovered upon expansion of MC@NLO results in a s . d) Hard emis¬ 
sions are treated as in NLO computations, e) Soft/collinear emissions are 
treated as in MC. f) The matching between hard- and soft-emission regions 
is smooth, g) MC hadronization models are adopted. Condition c) is the 
requirement that the MC@NLO be not affected by double counting prob¬ 
lems. In ME corrections, double counting corresponds to generating the 
same kinematical configuration twice (by hard process generation, and by 
means of the shower). On the other hand, the presence of negative weights 
in the MC@NLO implies that a double counting can be due both to an 
excess and to a lack of emission. 

The implementation of the MC@NLO formalism is technically compli¬ 
cated in QCD. However, it is possible to understand the basic features of the 
findings of ref. Q by studying a toy model. In this model, a system can ra¬ 
diate massless “photons”, whose energy I denote by x, with 0 < x < x s < 1, 
x s being the energy of the system before the radiation. After the radiation, 
the energy of the system is x' s = x s — x. The system can undergo several 
further emissions; on the other hand, one photon cannot split further. The 
NLO predictions for the mean value of an observable O (possibly including 
0 functions that define an histogram bin) is 
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where a is the coupling constant, and B, V, and R(x) are the Born, (fi¬ 
nite) virtual, and real contributions to the cross section respectively. The 
function O(x) returns the value of the observable O in the cases of real 
photon emission (x ^ 0), and of virtual photon emission or of no emission 
(x = 0). Eq. is based on standard subtraction techniques for the cancel¬ 
lation of infrared divergences. In the toy model, a MC can be defined with 
the following Sudakov form factor 


A(xi,x 2 ) = exp 
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where Q[z ) is a monotonic function, such that 0 < Q(z) < 1, and 
]im z —>o Q(z) = 1- The form factor in eq. dD can be used in a standard 
way to generate showers in which the system can emit any number of pho¬ 
tons. Since this number is in general larger than 1, the function 0(x ) is 
not sufficient to describe the observable O; we denote the distribution in the 
observable O obtained from the MC by 


Imc{0, x m (x)) , 


(3) 


where x M (x) indicates the initial condition for the shower (the system has 
energy 1 — x). 

A naive attempt at matching NLO computations and MC simulations 
amounts to formally replacing 0(x) with I MC (0,x M (x)) in eq. (Tj): 
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Here, a couple of MC runs are performed, one in which the initial condition 
is given by the system plus a photon of energy x x M (x))), and one 

in which there is no photon emission prior to the shower (I MC (0,1)); the 
results are then weighted with aR{x)/x and B + aV — aB/x respectively, 
as in eq. (0). Unfortunately, this does not work: the procedure in eq. (Q) 
does have double counting. The solution proposed in ref. (lj] (called modified 
subtraction there) is the following: 
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which is obtained from eq. by subtracting and adding the quantity 
Imc(0,x m ) aBQ(x)/x, using x M = x M (x) in the first and x M = 1 in the sec¬ 
ond term introduced in this way. These two terms compensate, at O(a), the 
shower contribution due to BI MC (0, 1), and therefore double counting does 
not occur. This happens because these terms coincide with the 0(a) term in 
the expansion of the Sudakov (3), the two different kinematics (x M (®) and 
1) accounting for the emission probability and for the no-branching proba¬ 
bility. Furthermore, the weights appearing in eq. (|5|) are finite (this is not 
the case in eq. ([§)), and this allows efficient unweighting. 

In order to deal with the case of QCD, one has to formally translate 
eq. (|5|). This is essentially a technical task; however, there are a couple 
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conceptual issues which do not arise in the simple toy model. In QCD, 
there are two types of non-UV singularities (soft and collinear), whereas 
in the toy model only soft singularities appear. In the case of collinear 
emission from initial-state partons, the resulting kinematics is one in which 
one massless parton is going down the beam line with non-zero energy. Such 
a configuration cannot be used as an initial condition for a parton shower; 
however, one can always boost the system to a frame in which the emitted 
parton has zero energy. MC-wise, this procedure can be made efficient by a 
suitable redefinition of the Bjorken cc’s |]]. The second problem is that, in 
the soft limit, the 0(a s ) term in the expansion of the MC result does not 
have the angular distribution expected from QCD factorization theorems, 
and this prevents a straightforward use of this result in the construction of 
the MC@NLO. However, the total amount of soft energy radiated by the 
MC is in agreement with QCD expectation, and this is sufficient for any 
infrared-safe observable to be predicted correctly. 

A public version of the MC@NLO advocated in ref. [1]] has been pre¬ 
sented in ref. [j|]. The production of standard model vector boson pairs ( ZZ , 
W + W~, and W±Z) in hadronic collisions has been implemented, matching 
the relevant NLO computations with Herwig. The code can be downloaded 
from tittp: //www.hep.phy. cam. ac .uk/theory/webber/MCatNLO. The im¬ 
plementation of the hadroproduction of heavy quarks in under way [tf]. 
Apart from its phenomenological relevance, it has to be stressed that the 
solution of the technical problems in the implementation of this process will 
make the implementation of any other production process in the MC@NLO 
a straightforward task. 

It is a pleasure to thank B. Webber for numerous discussions on this and 
other matters. 
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